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ABSTRACT 

We present the results of CO ( J = 3 — 2) and CO ( J = 1 — 0) mapping observations 
toward the active cluster forming clump, L1688, in the p Ophiuchi molecular cloud. 
From the CO ( J = 3 — 2) and CO ( J = 1 — 0) data cubes, we identify five outflows, 
whose driving sources are VLA 1623, EL 32, LFAM 26, EL 29, and IRS 44. Among the 
identified outflows, the most luminous outflow is the one from the prototypical Class 
source, VLA 1623. We also discover that the EL 32 outflow located in the Oph B2 
region has very extended blueshifted and redshifted lobes with wide opening angles. 
This outflow is most massive and have the largest momentum among the identified 
outflows in the CO ( J = 1 — 0) map. We estimate the total energy injection rate due to 
the molecular outflows identified by the present and previous studies to be about 0.2 L Q , 
larger than or at least comparable to the turbulence dissipation rate [~ (0.03 — O.1)L0]. 
Therefore, we conclude that the protostellar outflows are likely to play a significant role 
in replenishing the supersonic turbulence in this clump. 
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Introduction 



It is now well accepted that the majority of stars in our Galaxy form in clusters (jLada &; Lada 



2003). Recent observations have revealed that active cluster forming regions are not distributed 
uniformly in parent molecular clouds, but localized and embedded i n dense clumps with typical 
sizes of ~ 1 pc and masses of 10 2 ~ 1O 3 M0 (e.g., iRidge et al.ll2003l ). More than half of all the 
young stellar populations as sociated with a parent molecular cloud are confine d in such pc-scale, 



cluster forming clumps (e.g. jLada et al. 



1991 



Carpenter 



2000 



Allen et al.ll2006l ) , where many stars 



at different stages of formation are born in close proximity to one another. For example, in the 
Perseus molecular cloud, one of the nearby well-studied star forming regions, two pc-scale cluster 
formin g clumps, IC 348 and NGC 1333, contain about 80% of the young stars associated with this 
cloud ([Carpenter! 12000 ) . The mass fraction of molecular gas occupied by these two regions in the 
parent molecular cloud is less than a few tens %. In such dense clumps, the feedbac k from the older 
generation of stars affects the formation of the y ounger generation of stars (e.g., iNorman Silk 
1980l : lMcKeel[l989l : IVazauez-Semadeni et allboioh . 



Among the stellar feedback processes, protostellar outflows have been considered to be one 
of the impo rtant mechanism s to control the structure and dynamical properties of cluster form- 
ing clumps (jBallv et al.lll996l : Ishu et al.ll200ol : Matzner fc McKeelhoOol : lMatznerll2007r ) because the 
outflows from a group of young stars interact with a substantial volume of their parent clump by 
sweeping up the gas into shells. Indeed, recent numerical simulations of cluster formation have 
demonstrated that the protostellar outflo ws largely regulate the structure formation and star for 



2010 



mation in a dense c l uster forming clump (ILi &: Nakamurall2006l : iNakamura Lil 120071 ; IWang et al 



Li et al.ll2010i ) . ILi &: Nakamural (120061 ) showed that the supersonic turbulence in dense clumps 



can b e maintained by the momentum injection from the protostellar outflows (see also lCarroll et al 



20091 ). and thus the clumps as a whole can be supported by the turbulent pressure due to the proto- 
stellar outflows_jJgainst global gravitational collapse at least for several dynamical times. Further- 



Nakamura Lil (J2007J) showed that the global star formation efficiency tends to be reduced 



more, 

by the momentum injection from the protostellar outflows, although local star formation can often 
be triggered by the dynamical compression due to the protostellar outflows. These studies indicate 
that the protostellar outflow feedback is a key to understanding the formation of stars in a pc-scale 
cluster forming clump. 

Recent millimeter and submillimeter observations of nearby pc-scale cluster forming clumps 
have revealed the significant role of protostellar outflows in clustered star formation. An excellent 
example is NGC 1333, where several protostellar outflows apparently overlapping and interacting 
with themselves are found in the central region of the dense clump; those outflows infl uence the 



distribution of the dense gas significantly (jSandell Knee 



2001 



Gutermuth et al]l2009T ). The less 



dense parts, seen as cayities, are found to be filled with 



(Kn ee h Sandell 



2000 



Walsh et al 



2007 



Hatchell et al 



basis of the CO ( J = 1 — 0) observations, iQuillen et al 



high-velocity gas created by the outflows 



20071: lHatchell & Dunhamll2009h . On the 
found 22 cavities with sizes of 0.1 



0.2 pc in diameter in the central region of NGC 1333 (see also ICunningham et al.ll2006l ). Those 
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cavities appear to expand slowly. They interpreted those cavities as remnants of past YSO outflow 
activity, although the role of those "fossil" cavities in cloud dynamics is still unclear. 

Another example is the Serp ens clump, where several powerful outflows were identified by 
the CO (J = 2 — 1) observations (jDavis et al.lll999l ). The energy injection rate of these outflows 
is comparable to or somewhat larger than the dissipation rate of the turbulent energy in this 



clum p, indicating that the o utflows are the main source of the supersonic turbulence (|Sugitani et al 



20101 ) . ISugitani et al.l (|2010l ) suggested that the outflows in this clu mp have enough mo mentum to 
support the whole clump against global collapse. For NGC 2264-C, iMaurv et al.l (|2009l ) suggested 
that the parent clump cannot be sufficiently supported by the outflows in the current generation, 
although they are likely to contrib ute much to ma i ntaini ng the observed turbulence. In the Orion 
Molecular Cloud-2 FIR 3/4 region. IShimajiri et al.1 (|2008l ) found an observational example showing 
that the dynamical interaction of protostellar outflows with a dense clump may have triggered 
fragmentation of the clump, leading to future star formation. In this paper, to gain further evidence 
of the significant role of protostellar outflows in clustered star formation, we investigate the outflow 
activity in the dense part of the p Ophiuchi molecular cloud (hereafter the p Ophiuchi main cloud) , 
one of the nearby pc-scale cluster forming clumps, on the basis of CO (J = 3 — 2) and CO (J = 1 — 0) 
observations. 

The p Ophiuchi main cloud is a suitable object to investigate the star formation process in a 
cluster forming clump because this cloud is the nearest pc-scale cluster forming clump at a distance 
of about 125 pc (e.g., Lombardi et al. 2008; Loinard et al. 2008; Wilking et al. 2008). The p Oph 
clump is known to harbor a rich cluster of young stellar objects (YSOs) in different evolutionary 
stages (Wilking et al. 2008; Enoch et al. 2008; J0rgens en et al. 2008) . On t he basis of the recent 
infrared observations using the Spitzer space telescope, iPadgett et al.l (|2008l ) identified about 300 
YSOs, including Class 0/I/II/III sources, most of which are associated with the central region of 
the p Ophiuchi main cloud. Furthermore, recent millimeter and submillimeter dust continuum 
observations of p Oph have revealed that a large number of dense cores are concentrated in the 
central dense part of the clump and their mass spectra ar e similar in shape to the stellar IMF 
(jMotte et al.lll998l : iJohnstone et al.ll2000l : IStanke et al.ll2006l ). These facts indicate that active star 
formation is still ongoing in this cloud. 

So far, more than 10 molecular outflows from the YSOs have been identified in p Oph (e.g., 



Bontemps et al.lll996l ; ISekimoto et al.lll997l ; iKamazaki et al 



2003 



Wilking et al.ll200§). The most 



spectacular outflow in p Oph is the one from a prototypical Class source, VLA 1623, (jAndre et al 



1990 



Dent et al.lll995l ). which is associ ated with H e rbig-H aro objects, many knots of the H2 emis- 



sion at 2.12 /im, and the water masers. iDent et al.l (|1995l ) mapped the molecular outflow from the 
VLA 1623 using CO( J = 3 — 2) emission. They revealed that the outflow is highly -collimated with 



an op ening angle less than 1.6 degree, extending over a length of about 0.5 pc. IKamazaki et al 



(|2003l ) discovered three other molecular outflows in the Oph A and B2 regions, a lthough they did 
not su ccessfully map the entire outflow lobes. In the southern part of p Oph, iBussmann et al 



( 20071 ) mapped the outflows from EL 29 and LFAM 26. These outflows appear to be less powerful 
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than the VLA 1623 outflow. However, the physical properties of the molecular outflows previously 
identified have not been well determined because of the lack of the large-scale mapping observa- 
tions. To study how the outflows affect the clump dynamics and local star formation, we have 
performed the CO ( J = 3 — 2) and CO ( J = 1 — 0) observations toward the p Oph main cloud, 
using the ASTE 10 m and Nobeyama 45 m telescopes, respectively, and have derived the physical 
properties of the outflows identified in the clump. 

The rest of the paper is organized as follows. First, we present the details of our CO (J = 3 — 2) 
and CO ( J = 1 — 0) mapping observations in Section [2 In Section 3, we derive the physical 
parameters of the molecular outflows identified from our CO (J = 3 — 2) and CO (J = 1 — 0) 
maps. Then, we make a brief discussion of how the ob served outflow s inter act with the dense 



cores identified from the H 13 CO+ (J = 1-0) map by IMaruta et al.1 (|2010l ) in Section [U We 
also estimate how the observed outflows contribute to turbulence generation in this clump, and 
compare the clump with other nearby cluster forming clumps, Serpens and NGC 1333. Finally, we 
summarize our main conclusion in Section 



2. Observations 
2.1. CO ( J = 3 - 2) Observations 

The CO (J = 3 — 2; 345.7959 899GHz) mapping observations were carried out with the ASTE 



10 m telescope (jEzawa et al.ll2004l ) during the period of 2004 October. The beamsize in HPBW 
of the ASTE telescope is 22", which corresponds to 0.013 pc at a distance of 125 pc. The main- 
beam efficiency, 7732, was 0.35 at 345 GHz. We used a 345 GHz SIS heterodyne receiver, which 
had the typical system noise temperature of 175-350 K in DSB mode at the observed elevation. 
The temperature scale was determined by the chopper-wheel method, which provides us with the 
antenna temperature, T5, corrected for the atmospheric attenuation. As a backend, we used four 
sets of 1024 channel autocorrelators, providing us with a frequency resolution of 125 kHz that 
corresponds to 0.11 km s _1 at 345 GHz. The position switching technique was employed to cover 
the dense region of the p Ophiuchi main cloud, about 23' x 23' area, corresponding to about 0.8 
pc x 0.8 pc at a distance of 125 pc. After subtracting linear baselines, the data were convolved 
with a Gaussian-tapered Bessel function (Magnum et al. 2007) and were resampled onto a 5" grid. 
The resultant effective FWHM resolution is 40". The typical rms noise level is 0.28 K in brightness 
temperature (T m b) with a velocity resolution of 0.4 km s _1 . 



2.2. CO (J = 1 - 0) Observations 



As shown in the n ext section, our CO (J = 3 — 2) observations revealed that the outflow 
previously identified by iKamazaki et al.1 (120031 ) in the Oph B2 region, the EL 32 outflow, is very 
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extended. Unfortunately, the CO (J = 3 — 2) map does not cover the whole redshifted outflow lobe. 
To map the whole extent of the EL 32 outflow, we performed the CO ( J = 1 — 0) observations using 
the Nobeyama 45 m radio telescope. The CO ( J = 1 — 0; 115.2712018 GHz) mapping observations 
were carried out with the 25-element focal plane receiver BEARS equipped in the Nobeyama 45 
m telescope during the period from December 2009 to May 2010, as a part of the Nobeyama 45m 
Legacy Project (see http://www.nro.nao.ac.jp/). At 115 GHz, the telescope has a FWHM beam 
size of 15" and a main beam efficiency, r/io, of 0.32. At the back end, we used 25 sets of 1024 channel 
auto-correlators ( ACs) which have bandwidths of 32 MHz and frequency resolutions of 37.8 kHz 
(jSorai et al.ll2000l ). The frequency resolution corresponds to a velocity resolution of 0.1 km s -1 at 
115 GHz. During the observations, the system noise temperatures were in the range between 300 to 
900 K in DSB at the observed elevation. The standard chopper wheel method was used to convert 
the output signal into the antenna temperatures (TJ), corrected for the atmospheric attenu ation. 
Our mapping observations were made by the OTF mapping technique (jSawada et al.ll2008l ). We 
adopted a spheroidal function as a gridding convolution function (GCF) to calculate the intensity 
at each grid point of the final cube data with a spatial grid size of 12". The final effective resolution 
of the map is 30". The rms noise level of the final map is 1.0 K in brightness temperature (T m b) at 
a velocity resolution of 0.4 km s . 



3. Results 

3.1. CO ( J = 3 - 2) and CO ( J = 1 - 0) maps 

In Figures [T] and El we present the CO ( J = 3 — 2) and CO ( J = 1 — 0) integrated intensity 
maps in the ranges of ulsr = —9.8 ~ +12.2 km s _1 and t>LSR = —3.2 ~ +8.0 km s _1 , respectively. 
In Figures [3] and HI the distributions of the blueshifted and redshifted CO ( J = 3 — 2) and CO (J 



1 — 0) components, respectively, are also shown on the 850/mi images obtained by I Johnstone et al 



(|2000l ). In Figure EK, several dense subclumps, which are recognized in the 850 /mi images, are 
designated by A, Bl, B2, C, D, E, and F. By comparing the H 13 CO + spectra with the 12 CO 
(J = 3 — 2 and J = 1 — 0) ones at the same position, the shapes of the 12 CO (J = 3 — 2) and 
12 CO ( J = 1 — 0) line profiles appear to be caused mainly by self-absorption rather than multiple 
velocity components (see also Figures [5] and [6]) . 

The CO ( J = 1 — 0) map in Figure [2] shows a clear break between the east and west parts of the 
p Oph clump, although the radial velocities of the two parts are not so different (t>LSR ~ 3.5 — 4.5 
km s _1 ). The CO ( J = 3 — 2) integrated intensity map also appears to reasonably follow that of 
CO (J = 1 — 0), although the detailed distributions are different. The strongest CO ( J = 3 — 2) 
emission comes from the western side of Oph A (see Figure H]). As shown later, the region with the 
strongest CO ( J = 3 — 2) emission has both the blueshifted and redshifted components that are 
likely to be due to the outflow from the VLA 1623 source. In the entire Oph A filament, the CO 
( J = 3 — 2) component blueshifted from the systemic velocity of Oph A is predominantly strong, 
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suggesting that Oph A is compressed from the far side by the photodissociation region (PDR) 
excited by a young B3 star, SI. 

The second strongest CO (J = 3 — 2) emission comes near the north-west edge of the ob- 
servational box, in which the CO (J = 1 — 0) integrated intensity is the strongest. This region 
has both broadly distributed blueshifted and redshifted components, and may be influenced by the 
PDR excited by the B2 star, HD 147889. There are two other areas having moderately strong CO 
( J = 3 — 2) emission. One is located near the north-east part of the CO ( J = 3 — 2) observational 
box, and the other is an elongated structure running from the box center to the south-east part. As 
discussed below, these two components originate from the high velocity gas by molecular outflows. 



3.2. Outflow Identification 



Using the CO (J = 3 — 2) and C O ( J = 1 — 0) data cubes, we identified molecular outflows 



as follows (e.g., iTakahashi et al.ll2008l ). First, we scrutinized the velocity channel maps of both 
the data to find localized blueshifted or redshifted emission in the area within a radius of 1' from 



each Y SO. Here, we used the YSO list obtained by I Jorgensen et al.1 (120081 ) and Ivan Kempen et al 
(|2009i ). From our CO ( J = 3 — 2) data, we identified 4 outflows. All these outflows can also be 
clearly recognized in the CO ( J = 1 — 0) map. In addition, we identified one outflow from the CO 
( J = 1 — 0) data, which is located outside our CO ( J = 3 — 2) map. All the outflows identified in 
the present paper were already detected by previous studies. The characteristics of the identified 
outflows are summarized in Table [TJ In Figures [5] and [6l we present the CO ( J = 3 — 2) and CO 
( J = 1 — 0) spectra centered on some of the high-velocity lobes, respectively. Each of the spectra 
shows signs of significant self-absorption due to the ambient gas at ambient velocities; both the 
CO line emission is clearly optically thick. In the following, we discuss the characteristics of the 
outflows identified in the present paper, and then derive their physical parameters. 



3.3. Individual Outflows 



3.3.1. VLA 1623 



In Figure O two highly-collimated, strong blueshifted-components labeled (a) and (b) can be 
clearly recognized. The y trace the collimat ed (south-east) outflow from the VLA 1623 source, which 
was first discovered bv lAndre et al.l (|1990l ) with the CO ( J = 2 — 1) observations. The north-west 
lobe of this outflow, labeled (c), is less prominent than the south-east one. The north-west flow 
has the low-velocity blue and redshifted components, suggesting that the outflow axis is almost 
p arallel to the plane of the sky. This is in good agreement with the CO ( J = 2 — 1) observations 



of I Andre et al.1 (Il990l ). However, as shown in the next subsection, some physical quantities of the 



outflow are significantly different from theirs. This is due to the fact that their observations did 
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not cover the whole VLA 1623 outflow. For example, our outflow size is larger by a factor of 5, the 
outflow velocity is smaller by a factor of 2, and thus the estimated dyn amical time of the so uth-east 
component is longer by a factor of about 10 than that estimated by lAndre et al.l (|1990l ). When 
the outflow parameters are esti mated using the sam e area as they mapped, our quantities agree 
reasonably with those derived bv lAndre et al.l (|1990l ). 



The whole outflow lobes [the components (a), (b), and (c)] were first mapped bv lDent et al 

(|1995l ) using the CO (J = 3 — 2) observations. They found that the blueshifted flow is extended over 
a length of about 0.5 pc with a narrow opening an g le of less than 1.6 degree. Our CO ( J = 3 — 2) 



map is in good agreement with that of lDent et al.l (|1995l ). Since this component is distinct in our 
CO (J = 1 — 0) map presented in Figured! we think that it can be detected with the CO ( J = 2 — 1) 
emission . Our total mass and momentum of the outflow are reasonably in good agreement with 
those of Dent et al. (jl995l ). taking into account the different adopted fractional abundance ratio 
and distance to p Oph, which cause the difference in the mass estimation by a factor of 3. The 
blueshifted components are also remarkable in the line profiles presented in Figures [5] and [6j in 
which the line profiles of the components (a) and (b) show strong blueshifted wings. 



In Figure [7J, we sho w the positions of the H2 knots identified by iGomez et al.l (|2003l ) and 
Khanzadyan et al.1 (120041 ) . Many H2 knots are associated with the blueshifted lobe labeled (a) in 
Figure 03- These H2 knots are clearly recognized in the IRAC channel 2 (4.5 //m) image taken by 
the Spitzer space telescope (IZhang fc Wa ng 2010j)- iKhanzadyan et al.l (J2004J) compared their H2 
knots with those identified by Gomez et al. ( 20031 ) to measure the proper motions of the H2 knots, 
and estimated the flow speed of about 130 km s _1 , which is much faster than the outflow speed of 
about 20 km s _1 estimated from our CO ( J = 3 — 2) data. 



3.3.2. Elias 32 and IRS U 



In our CO ( J = 3 — 2) image shown in Figure O we can clearly recognize a spatially extended 
outflow in the Oph B2 region (labeled (d) and (e)). The line profiles presented in Figures [5] and [6] also 
show distinct high-velocity wing emission for both the blueshifted and redshifted components. In 
particular, the blueshifted wings of both the CO (J = 3 — 2) and CO (J = 1 — 0) line profiles are very 



strong . This outflow appears to be from a Class I YSO, Elias 32, as suggested by lKamazaki et al 
(|2003l ). Although in our CO (J = 3 — 2) image, most of the redshifted component labeled (e) 
appears to be outside the observed area, our CO (J= 1 — 0) map have successfully revealed that the 
redshifted flow also extends with a wide opening angle toward the south-east direction (see Figure 

U). 

The blueshifted (north-west) lob e appears to b e anti- correlated with the dense gas detected in 
H 13 CO + emission in the B2 region by lMaruta et al.l ( 2010 ) (see Figure EK). Around the peak of the 
CO ( J = 3 — 2) emission in the blueshifted lobe, the H 13 CO + emission is less prominent, indicating 
that the dense gas is swept out of the region. In addition, the dense gas well follows the south-east 
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part of the blueshifted lobe, showing an arc- like shape. According to iMaruta et al.l (|2010l ). a small 
hole is seen in their H 13 CO + data cube in the southern part of the Oph B2 region (labeled No. 
1 in their Figures 2b and 2c), which is located near the north-west edge of the redshifted lobe. 
Two other holes are also seen near the northern (No. 2 in their Figure 2b, corresponding to the 
south-east edge of the blueshifted lobe) and southern (No. 3 in their Figure 2b, corresponding to 
the north-west edge of the redshifted lobe) edges of the B2 region. These holes are likely to be 
created by the EL 32 outflow. Furthermore, t he blueshifted lobe is extended beyond the B2 region, 
reaching the No. 4 and 5 arcs in Figure 2b of IMaruta et al.l (|2010l ). These structure s suggest that 



the O ph B2 region is strongly influenced by the EL 32 outflow. See Section 3 of IMaruta et al 
(|2010l ) for more detail. 



Near the southern boundary of the observed area with CO ( J = 1 — 0) presented in Figure 21 
which is outside the CO ( J = 3 — 2) map, we can see a weak blueshifted and redshift ed lobes, which 



are p robably the outflow lobes from IRS 44. This outflow was first discovered by iTerebev et al 
(|1989l ) by means of the C O (J = 1 — 0) interferometric observations using the Owens Valley 



Millimeter Interferometer. iJorgensen et al.l (|2009l ) presented the three-color Spitzer IRAC image 
of several knots that are proposed to be associated with the outflow from IRS 43 or other nearby 
YSO including IRS 44 (see their Figure B.l). Although these knots are clearly along a line pointing 
directly back toward both IRS 43 and 44, the direction of our redshifted lobe of the IRS 44 outflow 
does not match that suggested by the knots. Thus, these knots are probably not associated with 
the IRS 44 outflow. 



3.3.3. Elias 29 and LFAM 26 



In the southern part of the CO ( J = 3 — 2) map (Figure [3]), we can find 4 redshifted (labeled 
(f), (g), (h), an d (i)) and 1 blueshifted (labeled (j)) components. This region was mapped by 



Bussmann et al. J 2007h using the C O (J = 3 — 2) emission, but our map covers a larger area 
than that of iBussmann et al.l (120071) . The over all feature of the CO (J = 3 — 2) map is in good 
agreement with that of IBussmann et al.l (|2007l ). although we found a new redshifted component 
labeled (i) in the southern part of this region. The redshifted components labeled (h) and (i) and 
the blu eshifted component lab eled (j) are most likely to be from a Class I YSO, EL 29. As pointed 
out by IBussmann et al.1 (120071). the E L 29 outflow axis has an S-like shape, traced by several H2 
knots detected by lGomez et al.1 (120031 ) and lYbarra et al.1 (120061 ). In addition, our CO ( J = 3 — 2) 
and CO ( J = 1 — 0) maps in Figures [3] and [4] show that the EL 29 outflow has both the blueshifted 
and redshifted components on the north side, suggesting that the outflow axis is almost parallel to 
the plane of the sky. The reason why the EL 29 outflow is so deformed is unclear. One possibility 
is the interaction with the magnetic fi eld. Recently, the magnetic field structure toward p Oph has 
been mapped by iTamura et al.l (|2010l ). They found that the ordered magnetic field penetrates the 
Oph E and F filaments almost perpendicularly. The direction of the magnetic field coincides well 
with the outflow axis of the redshifted component. Another possibility is the dynamical interaction 
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with dense gas. In fact, the redshifted outflow lobe is located in the less dense region between the 
Oph E and F regions. This implies that the outflow hit the dense gas ridge and the direction was 
altered. Both effects might be responsible for the shape of the redshifted outflow lobe. 

On the other hand, we interpret that the other redshifted components (f ) and (g) are associated 
wit h the outflow from a Class I YSO, LFAM 26, because a line connecting the two H2 knots detected 
by iKhanzadvan et al.l (|2004l ). f04-01a and f04-01b, reasonably coincides with a line between the 
outflow components (f) and (g). Furthermore, the bow-shocked-like shapes of the two knots support 
our interpretation that the component (g) is coming from LFAM 26. As seen in the CO ( J = 1 — 0) 
map in Figure [H the redshifted components labeled (f) and (g) are also associated with weak 
blueshifted components, suggesting that the LFAM 26 outflow axis is nearly parallel to the plane 
of the sky. This interpretation is consistent with the fact that the disk around LFAM 26 is almost 
edge-on with the rotation axis parallel to the east- west direction (jDuchene et al.1 12004 ) . Although 
the VLA 1623 outflow is still the most powerful in this region, the estimated speeds of the EL 29 
and LFAM 26 outflows are relatively as large as that of VLA 1623 when the inclination angles of 
the three outflows are taken into account. 



3.4. Derivation of Physical Quantities 

In the above subsection, we identified 5 outflows from the CO ( J = 3 — 2) and CO (J — 1 — 0) 
observations. Here, we derive the physical quantities of the 5 clear outflows, following the procedure 
described below. 

Under the assumption of local thermodynamical equilibrium (LTE) condition and optically 
thin emission, the molecular outflow masses derived from CO ( J = 3 — 2) and CO ( J = 1 — 0) can 
be calculated, respectively, as 

M 32 = M 32, 3 , (1) 
3 

and 

Mio = M io,j • (2) 

3 

Here, M^j and M\q j are the outflow masses at the j-th channel, derived from CO ( J = 3 — 2) 
and CO (J = 1 — 0), respectively, and are given by 

»(S) J (^)^. 



and 

Afioj = 3.4x10 



Vio- 4 y Vi25pcy 



arcsec 
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Vio 
0.32 



£^■(1-0) At, 
K km s _1 



Mr. 



J 



(4) 



where i denotes the grid index on the j-th channel. Here, the same excitation temperature T GX 
of 30 K is adopted for the both lines. The fractional abundance of CO relative to H2, Xqq, 
and the distance to the cloud, D 1 are adopted to be 1 x 10~ 4 (a typical value of high-density 



regions of molecular clouds; e.g., iGarden et al.l ()199ll )) and 125 pc (e.g., lLombardi et al 



2008 



Loinard et al.ll2008l : IWiiking et al.ll2008l ). respectively. The main beam efficiencies of the ASTE and 
Nobeyama 45 m telescopes, are r/32 = 0.35 and 7710 = 0.32, respectively. The integrated intensities 
of the CO (J = 3 - 2) and CO (J = 1 - 0) lines at the i-th grid are r/^ 1 £V T* Ai j (?> - 2)Av and 
^10 J2j — 0)Av, respectively, where we integrated the antenna temperatures above the 3 a 

noise levels in the velocity range given in Tables [2] and [3] (see Sections 12.11 and 12.21 for the values of 
the 1 a rms noise). The symbols AO and Av (= 0.4 km s _1 ) indicate the grid size and the velocity 
resolution of the data. 

For both the lines, the outflow mass M out , the outflow momentum -P ut> and kinetic energy 



E ou t are estimated as 



M 



(5) 



^Pout = ^2 M j\ Vj - V Bys I / cos £ , 



(6) 



and 



E, 



out 



y sys ) 2 /cos 2 e , 



(7) 



respectively, where Mj is the mass of the j-th channel, Vj is the LSR velocity of the j-th channel 
and V S y S is the systemic velocity of the driving source. We adopt V sys ~ 3.8, 4.3, 4.9, 4.9, and 4.2 
km s" 1 for VLA 1623, EL 32, EL29, LFAM 26, and IRS 44, respectively. Each systemic velocity 
was determined from the peak Vlsr of the H 13 CO + ( J = 1 — 0) emission averaged in the 2' x 2' 
area centered at the position of each source, where we used the H 13 CO + (J = 1 — 0) emission 
of the Nobeyama 45 m archive data whose fits data are available from the Nobeyama web page 



(http:/ /www.nro.nao.ac.jp/). The systemic velocitie s are in agreement wit h the values used by the 



previous studies: F svs ~ 3.7km s" 1 for VLA 1623 (IKamazaki et all 120031'). 3 9 km s" 1 for EL 32 
dKamazaki et al.ll2003l). 5 km s " 1 for EL 29 and LFAM 26 teussmann et al.ll2007h . and 3.8 km 
s _1 for IRS 44 ( Sekimoto et al.lll997l ). The angle £ is the inclination angle of an outflow, which 
is generally uncertain. The three outflows from VLA 1623, EL 29, and LFAM26 in the present 
paper are expecte d to be almost parallel to the plane of the sky. Therefore, we adopt £ = 80° (see 
Andre et al.lll99oh. For the remaining EL 32 and IRS 44 outflows, we adopt £ = 57.3°, following 
Bontemps et al.1 (|1996l ). The projected maximum size of each outflow lobe, .Robs; was measured at 
the 3 a contour in its channel maps. The maximum size of each outflow lobe was estimated by 
correcting for the projection with the inclination angle £ as i? ou t = Robs/ sin£. The velocity range of 
each outflow was determined using the channel maps: isolated components from high-velocity CO 
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emission are clearly identifiable at more than 3 a level within the velocity range, without significant 
contamination of emission from the ambient gas. We also note that for all the previously-identified 
outflows, the velocity ranges for the integration are more or less similar to those listed in Tables [2] 
and[3l 

From these quantities, we calculated the characteristic velocity V out = P ut/M out , dynamical 
time t dyn = R oxlt /V out , mass loss rate M out = M out /t dyn , outflow luminosity L out = E out /t dyn , 
and outflow momentum flux F ou t = -Pout/^dyn- Tables [2] and [3] summarize the outflow parameters 
derived from the CO (J = 3 — 2) and CO ( J = 1 — 0) data, respectively. The CO lines often become 
optically thick even toward the outflow wing components. Therefore, the physical parameters such 
as the outflow mass and momentum listed in Tables [2] and [31 where both the emission lines are 
assumed to be optically thin, give the lower limits. 

Most of the physical quantities estimated from CO (J = 3 — 2) and CO (J = 1 — 0) are 
in reasonably good agreement with each other. The total outflow luminosity estimated from CO 
( J = 3 — 2) is only about 6 % larger than that of CO ( J = 1 — 0). However, for the EL 32 outflow, 
the mass estimated from CO ( J = 1 — 0) is larger by a factor of 4 than that of CO (J = 3 — 2). 
This deviation causes large discrepancy in the other quantities such as the momentum, momentum 
flux, and luminosity. This deviation on the estimated outflow mass is greatest among those of the 
five outflows. This may imply that the EL 32 outflow is more evolved and thus the outflow gas 
density may be smaller than the critical density of the CO ( J = 3 — 2) transition, 3 x 10 4 cm~ 3 . 
If this is the case, we are probably underestimate the CO ( J = 3 — 2) outflow mass because of the 
LTE assumption. This interpretation is consistent with the fact that the dynamical time of the EL 
32 outflow is long and that the outflow lobes are most spread out with the widest opening angles. 



4. Role of Protostellar Outflows in Clustered Star Formation 
4.1. Interaction between the Outflows and Dense Gas in p Oph 

Recent observations suggest that the dynamical interaction of dens e gas with prot o stellar 



outflows is likely to trigger the formation of dense cores. For example, ISandell &: Kneel (|200ll ) 
performed the 850 pm dust continuum emission observations toward a nearby pc-scale cluster 
forming clump, NGC 1333, and found that the dense cores are distributed al ong multiple she l ls and 



filaments that are most likely created by protostellar outflows. Recently, IShimajiri et al.1 ((2008) 
observed the Orion Molecular Cloud-2 FIR 3/4 regions using the Nobeyama Millimeter Array, 
and revealed that the 12 CO emission toward the FIR 4 region shows an L-shaped structure in the 
position-velocity diagram, suggesting that an outflow from the FIR 3 region compressed the gas 
in the FIR 4 region, triggering the formation of multiple dense cores in the past few xlO 4 yr. 
These observations suggest that dynamical compression due to protostellar outflows regulates the 
formation of dense cores in pc-scale cluster forming clumps. Here, we discuss how the outflows in 
p Oph affect the structures of the dense gas. 
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Using the H 13 CO + ( J = 1 — 0) data obtained with the Nobeyama 45 m telescope, iMaruta et al 



(|2010l ) identified 68 dense cores in the central dense region of the p Ophiuchi main cloud. The 
H 13 CO + ( J = 1 — 0) transition has a critical density of about 10 5 cm -3 , and therefore is suitable 
as a dense gas tracer. In the following, we compare the H 13 CO + data with the outflows identified 
in the present paper, to examine whether the outflows are interacting with the dense gas in this 
region. 

The interaction between the CO outflows and the H 13 CO + dense gas seems distinct in the Oph 
B2 region. In Figure [8^,, we present the CO ( J = 3 — 2) images on the H 13 CO+ total integrated 
intensity map. We also show the positions of the H 13 CO + cores identified by IMaruta et all ( 201y) 
in Figure [8b by open squares. About 30 out of 50 H 13 CO + cores located in our CO ( J = 3 — 2) 
observed area are distributed around the outflow lobes. For example, in the Oph B2 region, the 
dense region detected by the H 13 CO + emission has an arc-like shape, following the south-east part 
of the blueshifted lobe. The part with strong CO (J = 3 — 2) emission in the blueshifted lobe 
appears to be less dense: the peak of the CO ( J = 3 — 2) blueshifted emission does not correspond 
to any peak of the H 13 CO + emission. 

The EL 32 outflow possibly injected some energy into the H 13 CO + dense cores. In Figure O 
we present t he line- widths o f the H 13 CO + cores against their radii, which is essentially the same as 
Figure 9a of IMaruta et alJ (12O10h . The H 13 CO+ cores associated with the Oph B2 region, which 
are indicated by the filled circles, t end to have larger v elocity widths. This tendency seems to be 
consistent with the observations by lFriesen et all (J2009J), who performed NH3 observations toward 
the Bl, B2, C, and F regions and found that the NH3 cores in the B2 region have somewhat larger 
velocity widths. In addition, they found that the Oph B2 region is slightly warmer (T ~ 15 K) 
than the other regions (T ~ 12 K). iFriesen et all (|2010l ) also revealed that the ratio of the NH3 
abundance to N2H + abundance in the Oph B2 region appears to be smaller than that of the other 
regions. These facts appear to support the idea that the Oph B2 region is largely affected by the 
outflow activity. However, there is a possibility that the cores in the B2 region simply tend to 
follow the line width-radius relation for which the cores with larger radii have larger line widths, 
althoug h the line-widths of the H 13 CO + cores seem to be almost independent of the core radii for 
p Oph (jMaruta et al.ll2010l ). We need further observational evidence to quantify the interaction of 
the dense gas with the outflow. 

In addition to the EL 32 outflow, the VLA 1623 outflow also seems to interact with the dense 
gas. The head of the collimated blueshifted lobe of the VLA 1623 outflow is located just between 
the Oph Bl and C regions, w hich might be broken up into the two parts by the outflow. The No. 
33 H 13 CO + core identified by IMaruta et all (|2010l ) appears to collide with the VLA 1623 outflow 
near the edge of the Oph C region. There are also about 10 H 13 CO+ cores that are distributed 
around the collimated blueshifted lobe of the VLA 1623 outflow. 



The redshifted lobe of the EL 29 outflow has an S-shape axis, going through the less dense area 
between the Oph E and F regions. This suggests that the outflow was presumably bent into an 
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S-shape as a result of the dynamical interaction with the dense filamentary ridge. The redshifted 
lo be has the two peaks labeled (h) and (i), near which two H 13 CO + cores (No. 34 and 39 in Table 1 
of lMaruta et al.l (|2010l )) are located. The difference between the radial velocity of the cores (~ 4 — 5 
km s _1 ) and velocity of the redshifted lobe ranges from 2 and 3 km s -1 . Thus, the external pressure 
exerted on the core surfaces is expected to be significant. In fact, the cores have the virial ratios of 3 
— 6, suggesting that they are not gravitationally confin ed but compressed b y the external pressure. 
Furthermore, the detailed virial analysis performed by iMaruta et al.l (|2010l ) indicates that most of 
the cores in p Oph seem not to be confined by their self-gravity, but to be confined by the ambient 
turbulent pressures. This trend is in good agreement with recent numerical simulations of cluster 
formation by Nakamura & Li (2010, in preparation). These characteristics appear to be consistent 
with the idea that the outflows trigger the formation of dense cores in the central part of p Oph. 
However, the core distribution projected on the plane of the sky may be just a coincidence. In the 
3D space, the cores may not interact with the outflow lobes. To directly verify the interaction with 
the outflows, we need additional observations such as using SiO and CH3OH emission, known as 
shock tracers. 



The importance of the external pressures on core dynamical sta tes has been sugg e sted even in 



the ne arby low-mass star forming regions such as the Pipe nebula (jLada et al.l 1200a ). lLada et al 



( 20081 ) performed a detailed virial analysis of the cores identified from the extinction map of this 
region, and found that all the cores have virial ratios larger than unity, suggesting that most of the 
cores are pressure-confined. However, it is unlikely that the large external pressur e is due to the 



proto stellar outflows because star formation activity is quite low in the Pipe nebula (jForbrich et al 



2009). Furthermore, the typical external pressure of the Pipe nebula (P/k B ~ 8 x 10 4 K cm 3 ) is 
nearly an order of ma gnitude sma ll er tha n that of p Oph (P/ks ~ 3 x 10 6 K cm -3 ), where fee is the 
Boltzmann constant. lLada et al.l (J2008J) proposed that the source of the external pressure is due 
to the self-gravity of the Pipe nebula itself. In any case, the previous and our studies indicate that 
the ambient turbulent pressures of the molecular clouds play a significant role in core dynamics for 
both clustered and distributed star formation. 



4.2. A Dynamical Role of the Outflows in p Oph 



4-2.1. Turbulent Dissipation and Generation 



On the basis of numerical simulations, iMac Low! (Il999l ) derived the approximated formula for 
the energy dissipation rate of supersonic turbulence as 



L 



turb 



1/2MAV 2 



(8) 



where /(~ 0.33) is the nondimensional coefficient, AV is the I D FWHM veloci ty width, and 
is the driving scale of supersonic turbulence (see equation [7] of lMac Low! ([1999)) . Applyin g this 
formula, we evaluate the rate of turbulent energy dissipation in the p Ophiuchi main cloud. lLoren 
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(|1989al ) estimated the total gas mass of the whole p Ophiuchi main cloud (corresponding to the R21, 
R22, R24, R25, and R26 regions in his Table 1) to be 883 M using the 13 CO (J = 1 - 0) emission, 
where the mass is recalculated by using our assumed distance of 125 pc to the cloud. The ID 
FWHM velocity width and the mean d iameter of p O ph (L1688) are estimated to be AV^ ~ 1.5 km 
s _1 and about 1.6 pc, respectively (see iLorenl 1 1989al ) . The rate of the turbulent energy dissipation 
is estimated to be 

Lturb - O.O3L , (9) 

where we took the cloud diameter as the driving scale of the turbulence. We note that the driv- 
ing scale of the outflow-driven turbulence is rather uncertain because the driving scale may vary 
significantly between outflows. In fact, the size of the outflows detected from our observations 



range from 0.1 pc to 0.6 pc. For example, iMaurv et al.l (|2009l ) ado pted the cloud di ameter as the 



driving scale of the outflow-driven turbulence. On the other hand, lArce et al.l (|2010l ) adopted the 
driving scale of 0.2 pc for NGC 1333, although they mentioned that this may be a lower limit 
because the most outflow lob es they identified in NGC 1333 have larger sizes. From a theoretical 
consideration. IMatznerl (|2007l ) derived the driving scale of 0.4 pc fo r the out fl ow-dr iven turbulence. 



Nakamura &: Lil (|2007l ) obtained a driving scale similar to that of IMatznerl (|2007l ) on the basis of 
the numerical simulations of cluster formation. If we adopt = 0.4 pc, then the energy dissipation 
rate increases up to Lturb — O.12L . 

From the analysis in the previous section, we can evaluate the total energy injection rate 
(mechanical luminosity) due to the protostellar outflows in L1688 as 



L 



tot 



0.2Lrr 







(10) 



where E OVLt) k and idyn.fc are the outflow energy and dynamical time of the fc-th outflow, respec- 
tively, and we used the quantities obtained from the CO ( J = 3 — 2) observation s. Here, we also 



included the contribution f r om small outflows det e cted in previous ob servations (jBontemps et al 



1996 



Sekimoto et al.l Il997l ; iKamazaki et al.l 120031 ; IStanke et al .1 120061 ) . The physical parameters 



of such outflows are listed in Table 01 where all the quantities have been scaled such that they 
correspond to an assumed distance of 125 pc. The estimated L tot is sensitive to the inclination 
angles assumed. If we adopt £ = 57.3° (random orientation) for all the outflows, the above value 
is reduced to L tot ~ 0.03L Q . The estimated L tot gives a lower limit of the outflow energy injection 
rate because we assumed the optically thin for the CO line emission and we neglected the contribu- 
tion from undetected faint outflows. Therefore, we conclude that the total outflow energy injection 
rate is larger than or at least comparable to the turbulent dissipation rate. 



4-2.2. Force Exerted by the Outflows 



To clarify how the outflows affec t the d ynamical state of the cloud, we assess the force balance 
in the cloud, following IMaurv et al.l (|2009l ). To prevent the global gravitational contraction, the 
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following pressure gradient is needed to achieve the hydrostatic equilibrium: 



dP< 



grav 



dr 



_ G M(rMr) 



(11) 



where M(r) is the mass contained within the radius r and we assume that the cloud is spherical. 
The effect of magnetic field is taken into account by the factor (1 — a~ 2 ) and a is the mass-to- 
magnetic flux ratio normalized to the critical value (e.g., Nakano 1998). Assuming the density 
profile of p oc r~ 2 , the pressure needed to support the cloud against the gravity is estimated to be 

GM(R) 2 



P no 

1 grav 



1 



a 



87TP 4 

The force needed to balance the gravitational force is thus evaluated to be 

F grav ~ ^R 2 P giav (R) = (1 - a- 2 ) . 



grav\ 



2R 2 



(12) 



(13) 



Adopting M(R) = 883M , R = 0.8 pc (|Lorenlll989al ). and a = 1.4, F g rav can be estimated to be 



grav 



1.3 x KT 3 M km s^yr -1 



(14) 



where we used the same a as that estimated by ISugitani et al.l (|2010i ) for the Serpens main core, 
a nearby cluster forming clump similar to p Oph, although this value is highly uncertain. The 
moderately strong magnetic field of a = 1.4 can reduce the gravitational force by a factor of 2. On 
the other hand, the total force exerted by the outflows in this region is estimated as 

^"out,fc 



Fi 



tot 



^dyn,fc 



0.12 x 10" 



5 M km s _1 yr _1 



(15) 



The force due to the outflows, F tot is one order of magnitude smaller than that needed to stop the 
global gravitational collapse, Fgrav Therefore, the force exerted by the outflows seems not to be 
sufficient for directly supporting the whole cloud against its gravitational collapse. However, we 
assumed that the emission lines are optically thin in deriving the outflow parameters, and therefore 
the estimated total force of the out flows gives the lower l imit. If we adopt the optical depth of 

out ~ u.i a iu 1VJ.Q a yr _1 , 

grav • If this is the case, the force exerted by the outflows may significantly 
contribute to the cloud support. The accurate estimate of the optical depth is necessary to clarify 
the dynamical role of the outflows on the cloud support. 



r - 5 for the CO (J = 3 - 2) line (jKamazaki et al J 12003 ) . F out ~ 0.7 x 10" 3 M ra km s" 1 
which is about a half of F t 



4-2.3. Comparison with other nearby cluster forming clumps 

Recent observations of nearby cluster forming regions have revealed the detailed structure of 
pc-scale dense clumps and star formation activity including the outflows. Here, using the same 
method as in Section T4.2.21 we compare the dynamical states of two nearby cluster forming clumps, 
NGC 1333 and Serpens, with p Oph. 
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Table [5] shows the global p roperties of the nearb y cluster forming clumps. For the Serpens 
clu mp, we adop t ed th e results of ISugitani et al.l (|2010l ) who used the outflow parameters estimated 
by iDavis et al.l (|1999l ). For NGC 1333 , we estimated the physical quantitie s using the outflow 
parameters obtained bv lSandell fc Kned (|200ll ) and the clump mass obtained bv lRidge et al.l (|2003l ). 
The outflow parameters listed in Table [5] are estimated under t he assumpti o n of optically thin. 
Therefore, they give the lower limits. Recently, for NGC 1333, lArce et al.l (|2010l ) performed a 
detailed analysis on the effect of outflows in the parent clump using the 12 CO (J = 1 — 0) and 13 CO 
( J = 1 — 0) lines, taking into account the effect of opacity. The virial mass is estimated from the 
total mass and mean FWHM velocity width as 



M V iR = 210a 



-i 



pc 



AV 

km s - 



Mr. 







(16) 



where the dimensionless parameter a is set to 5/3, corresponding to the density profile of p oc r~ 2 . 
The cloud masses, velocity widths, and cloud radii are adopted from the references indicated in the 
caption of Table For all the clumps, the energy injection rate of the outflows seems larger than 
or comparable to the turbulent dissipation rate. Thus, we expect that the protostellar outflows 
can power supersonic turbulence in these three pc-scale cluster forming clumps. Our analysis also 
indicates that for Serpens and NGC 1333, the outflows can exert the significant force against global 
gravitational collapse. However, for p Oph, the total force of the observed outflows appears not to 
be enough to support the whole cloud even if we correct the effect of opacity by a factor of ~ 5. 
In fact, the virial ratio is quite small, and the clump should collapse globally. We note that the 
estimated virial ratio of p Oph is insensitive to the adopted line. For example, if we use the C 18 
(J = 1 — 0) emission, the virial ratio is estimated to be about 0.2, where we adopted the cloud 
mass of 494 M^. the velocity width of about 1.6 km s , and the radius of about 0.3 pc from 
Tachihara et al.l (hood ). 



Why is the estimated virial ratio of p Oph so small? One possibility is that the current star 
formation go t less active t han b efore, and therefore the cloud turbulence dissipated significantly. 
According to I Evans et al.1 (|2009l ). the current star formati on rate appea rs to be smaller than the 
value averaged over the cloud lifetime for p Oph. In fact, iLorenl (|1989bl ) obtained a very shallow 
linewidth-radius relation for the p Oph clump, based on the 13 CO (J = 1 — 0) observations. 
The shallow linewidth-radius relation might suggest that the large-scale turbulent motions have 
dissipated significantly because of the lack of the large-scale energy injection events. If this is the 
case, the clump as a whole might be in the course of the large-scale gravitational collapse. Another 
possibility is that the p Oph clump as a whole has been compressed by a large-scale shock, which 
has enhanced the cloud gravitational energy, reducing the virial ratio of the clump significantly. 
Several previous studies suggest that the p Oph cloud has been compressed by an expanding shell 
created by stellar winds and supernova ex plosions in the Sco OB2 association, and star formation 



in this cloud has been externally t riggered d Vrbal 1 1 9771: iMeyers et al 



de Geus 


1992: 


Motte et al. 


1998: 



Tachihara et al 



200 



3T 



1985 



Loren fc Wootten 



1986 



The cloud morphology and kinematics 
appear to support this scenario. For example, the south-western side of the cloud has a sharp edge, 
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whereas the north-eastern side has a less dense elongated structure. This is consistent with the 
scenario that the shell passing through the cloud from the south-western side to the north-eastern 
side woul d compress t he clo ud. The expanding shell is also observed with the HI emission. In 
addition, iMotte et al.l (|1998l ) revealed that the dense cores as well as several YSOs appear to be 
distributed along a linear chain that is roughly perpendicular to the direction of shock propagation. 
They proposed that the large-scale shock has c ompressed the cloud to create a dense ridge that 
has fragmented into many cores to form stars. Meyers et al.1 ()1985l ) suggested that the velocity 
difference of < 10 km s _1 between optical absorption lines of CH and CH + is caused by the shock 
transmitted into the dense clump. This yields the cloud-crushing time of t cc ~ a few 10 5 yr, which 
is about 5 — 10 times shorter than the typical age of the YSOs in this clump. Here, the cloud 
crushing time t cc is the timescale for the cloud to be crushed by the shock transmitted into the 
cloud. According to numerical simulations of the interaction of a large-scale shock with a cloud, 
the shock moti ons are converted into t he cloud turbulence sufficiently within several cloud crushing 
times (see e.g.. iNakamura et al.ll2006l ). The large-scale shocks ar e sometimes considered t o be one 
of the important driving sources of the cloud turbulence (e.g., iMcKee &; Ostrikerl 120071 ) . For p 
Oph, the observed turbulent motions seem too weak to support the whole cloud against the global 
gravitational collapse, although the shock due to the large-scale expanding shell has passed over 
some ~10 t cc ago. The large-scale shock seems not to have been sufficient to supply the supersonic 
cloud turbulence in this clump. 



Using the 13 CO (J = 1 - 0) data, IPadoan et al.l (120091 ) derived the power spectrum of the 
turbulent field in NGC 1333 and claimed that there is no evidence of outflow-driven turbulence in 
NGC 1333. This result apparently contradicts our conclusion. However, the 13 CO emission can 
trace only tenuous inter-core gas with densities as low as 10 3 cm~ 3 . Thus, their power spectrum is 
likely to be contaminated by the ambient components that are not associated with the dense cluster 
forming clump whose typical density is 10 4 cm -3 . In addition, it is difficult to detect protostellar 
outflows from the 13 CO ( J = 1 — 0) line because of its weak emission at the high-velocity wings. 
The line also tends to be strongly optically thick toward cluster forming clumps. Therefore, the 
13 CO emission may not be suitable for searching the influence by protostellar outflows in molecular 
clouds. More careful analysis will be needed. Denser molecular gas tracers such as the C 18 Q 
(J = 1 — 0) line are likely to be more appropriate for this kind of analysis. In fact, iBrunt et al 



(J2009J) investigated the turbulent field in NGC 1333 on the basis of the PCA analysis using the 
12 CO, 13 CO, and C 18 emission, and found that the power index of the turbulent spectrum derived 
from C 18 is significantly shallower than those of 12 CO and 13 CO □, implying that the outflows 
are likely to be responsible for the turbulence generation in the pc-scale cluster forming clump (see 



Recently, I Carroll et al. performed 3D MHD simulations of protostellar outflow-driven turbulence and 

demonstrated that the measurements of the turbulence power spectrum using the PCA significantly underestimate 
the contribution of protostellar outflows to the power spectrum, preventing the detection of the characteristic driving 
scale of the outflow-driven turbulence. Therefore, a role of the outflow-driven turbulence on the larger-scale turbulence 
in molecular clouds should be clarified further in future. 
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also lArce et al.ll2010l ). This result supports our conclusion. 



5. Summary 

We summarize the primary results of the present paper as follows. 

1. We have performed the CO ( J = 3 — 2) and CO ( J = 1 — 0) observations toward the nearest 
pc-scale cluster forming clump, p Ophiuchi main cloud, and identified 5 outflows, whose driving 
sources are VLA 1623, EL 32, LFAM26, EL29, and IRS 44. We found that the EL 32 and EL 29 
outflows identified by the previous studies are larger than previously mapped. 

2. We estimated the physical quantities of the outflows. The most luminous outflow was found 
to be the one from VLA 1623. We also discovered that the EL 32 outflow has the largest mass and 
momentum among the identified outflows. The physical quantities estimated from the two lines 
reasonably agree with each other. 



3. We compared our outflow images with the H 13 CO + map obtained by lMaruta et al.1 ( 20101 ) . 
and found that 30 out of 50 H 13 CO + cores appear to overlap with the outflow lobes on the plane 
of the sky. This might imply that the outflows influence the dense cores significantly. 

4. We estimated the physical parameters of the identified outflows and measured the total 
outflow energy injection rate and momentum flux. The total outflow energy injection rate is 
significantly larger than the dissipation rate of the supersonic turbulence in this region. We conclude 
that the outflows can power the supersonic turbulence in this region. 

In contrast, the outflows do not have enough momentum flux to support the whole clump 
against the global gravitational contraction. However, this result is sensitive to the assumed optical 
depth for the CO lines. If the CO lines are optically thick, then the contribution of the outflow 
force to the cloud support is expected to be significant. The accurate estimate of the optical depth 
for the CO lines toward the outflows is needed for further clarification of the role of the outflows 
on the cloud support. 

5. Using the literatures, we applied the same analysis to the two nearby pc-scale cluster forming 
clumps: Serpens and NGC 1333, and obtained the total energy injection rates and momentum 
fluxes. Our analysis indicates that for both the cluster forming clumps, the outflows are likely to 
play a significant role in the global cloud support and turbulence generation. 

6. The virial ratio of p Oph is estimated to be only 0.2 even if the contribution of the moderately 
strong magnetic support against global collapse is taken into account. Such a small virial ratio might 
suggest that the current star formation is inactive, and thus energy injection from the forming 
stars is not sufficient to support the clump. Alternatively, the small virial ratio might be a result 
of the large-scale shock compression. The observed weak cloud turbulence may suggest that the 
large-scale shock due to supernovae cannot supply the supersonic turbulence significantly to this 
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cloud. Instead, the large-scale shock may have triggered further star formation by increasing the 
gravitational energy. 

This work is supported in part by a Grant-in- Aid for Scientific Research of Japan (20540228, 
22340040). 



A. Mass of Molecular Outflows 



In this Appendix, we derive an expression for the outflow mass calculated from the CO (J = 
3 — 2) and CO (J = 1 — 0) emission under the assumption of local thermodynamic equilibrium 
(LTE) with the populations of all levels characterized by a single excitation temperature, T ex . The 
CO (J = 3 — 2) and CO (J = 1 — 0) emission associated with the molecular outflows is also assumed 
to be optically thin. In this sense, our estimated outflow masses are the lower limit. 

For molecular gas in LTE, the optical depth for the transition from upper level, J + l, to lower 
level, J, is given by 



T v = 



cM J+ i,j 2 J + 3 

8lTV 2 2 J + 1 



1 — exp 



hv 



kBT C : 



njL 
dv 



(Al) 



where c is the speed of light, A/+i,J is Einstein's A coefficient for the transition from the upper 
level to the lower level, v is the frequency of the transition, ks is the Boltzmann constant, nj is 
the density of the molecule in the lower level, J, and L is the thickness along the line of sight 
[Aio = 7.203 x 10~ 8 s- 1 and A 32 = 2.497 x lO" 6 s^ 1 for CO (J = 1 - 0) and CO (J = 3 - 2), 
respectively] . 



The column density of the molecule in the lower level J is then estimated as 



Nj 



I 



njLdv 



8irv 3 2 J + 1 



where the relation dv 



c?Aj +1 ,j2J + 3 
v/c)dv is used. 



exp 



hv 
kBT e: 



J r v dv , 



(A2) 



The column density of the molecule in the lower level J is also expressed using the total column 
density, JV to t, as 



and 



Z(T c: 



2 J + 1 

Nj = mD Ntotexp 



J^(2J + l)exp 
j=o 

kBT cx 



hBJ{J + 1) 
kBT cx 



hBJ(J + l) 



hB 



1 hB 
1 + -t-^t- + 



kBT ex 
1 



3k B T ex 15\k B T ( 



hB 



(A3) 

(A4) 
(A5) 
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where Z is the partition function and B is the rotational constant (B = 57.8975 GHz for CO). 
From eq. jA2j . we thus obtain the expression for the total column density as 



Ni 



8tti/ 



tot 



c 3 A J+hJ (2J + 3) 



Z(T ex ) exp 



hBJ(J + l) 

&bT cx 



1 — exp 



hv 

kBT c: 



T v dv . (A6) 



The optical depth t u can be obtained from the expression for the antenna temperature at the 
frequency v: 

T* A = ^ [/(Tex) - /(T bg )] [1 - exp (-t u )\ « ^ [/(T ex ) - /(T bg )] r„ , (A7) 



and 



i 



(A8) 



exp(hu/kBT) — 1 

where r\ is the beam efficiency and T bg is the temperature of the background radiation. The beam 
dilution factor is set to unity. We adopt the cosmic microwave background radiation with T bg = 2.7 
K as the background radiation. Then, the value of /(T ex ) is much larger than /(T bg ) when T ex ~ 30 
K, the adopted value in the present paper. 

Therefore, we neglect the second term in the right hand side of eq. |A7j . /(T bg ). The column 
density at the i-th grid on the j-th channel is then given using the relation v = 2B(J + 1) by 

87rk B v 2 „ trri \hB(J+ l)(J + 2) 



tat, 



-Z(T ox )exp 



2k B T c: 



(A9) 



h^Aj +1 j(2J + 3) 

Thus, the outflow mass calculated from the CO ( J = 3 — 2) and CO ( J = 1 — 0) emission is given 
respectively, by 

M 32 = ^ M 32J 



and 



Mio = ^Mioj , 



(A10) 
(All) 



where 



and 



M-. 



32,j 



fim H X C Q^lD 2 iVtot.ij 



9.4 x 10- 10 f ^ 
VlO" 4 



J (l25pc 



xT ex exp 



16.6 K 



M 



io,j 



1 + 



1 /^ bR 



3 vr e: 



+ 



A8 

arcsec 

1 /T b „ 



15 \T ( 



0.35 

2 



1 fJ2i^A,i,j 



Av 



K km s 



A/ 



i 



9.2 x 10- 9 ( ^ 
xT ex exp 



~ L ' D \ 2 ( A9 \ 2 



VlO" 4 J 
5.53 K 



125pc/ \ arcsec J \0 



mo 



.32 



K km s _1 



1 (T h 



3 V T, 



1 /T v 



15 V T f 



Mr. 







(A12) 



(A13) 
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Here, i is the grid index on the j-th channel, 7732 and 7710 are the main beam efficiencies of the 
ASTE and Nobeyama 45 m telescopes, respectively, tbh 2 is the mass of a hydrogen molecule, Xqo 
is the fractional abundance of CO relative to H2, and /i is the mean molecular weight of the gas 
and is set to 2.33. Q is the solid angle of the object and D is the distance to the object. 
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Table 1. Molecular Outflows Identified in the Present Paper 



Name 1 


Region 


Reference b 


Characteristics 


VLA 1623 


Oph A 


1,2,3,6,7 


Class 0, highly-collimated, strong blueshiftcd lobe 


EL32 


Oph B2 


1,7 


Class I, very extended lobes, most massive and most luminous 


EL29 


Oph E 


1,3,4,5,8 


Class I 


LFAM 26 


Oph E 


1,4 


Class I 


IRS 44 


Oph F 


1,4,8 


Class I 



"Names of driving sources 

b ( 1) This work; (2 ) lAndre et all lll990l); (31 iBontemps et al.l lll996h; (4) iBussmann et al.l i r2007f) ; (5) ICeccarelli et al.l hOO'j) ; 
(6) iDent et al.l il995r i; f7) iKamazaki et ail l l2003t) ; f8) ISekimoto et al.l il997fi 



Table 2. Molecular Outflow Parameters Obtained From CO ( J = 3 - 2) 



Name 


Velocity Range 


Vout 


Pout 


M out 


-Pout 


^dyn 


Mout 


Pout 


Pout 


Position 




(km s- 1 ) 


(km s~ 1 ) 


(PC) 


(1O- 2 M ) 


(M km s _1 ) 


(10 4 yr) 


(Meyr- 1 ) 


(M km s _1 y r_1 ) 


(io- 2 p ) 




Blue Lobe 


VLA 1623 


-5.0 ~ 1.0 


21.7 


0.51 


2.2 


0.49 


2.3 


9.8 x 10" 7 


2.1 x 10" 5 


4.2 


S-E 


VLA 1623 


-5.0 ~ 1.0 


20.4 


0.098 


0.55 


0.11 


0.47 


1.2 x 10~ 7 


2.4 x 10~ 5 


4.3 


N-W 


EL 32 


-5.0 ~ 0.5 


9.73 


0.25 


2.3 


0.22 


2.6 


8.8 x 10~ 7 


8.6 x 10~ 6 


0.71 


N-W 


LFAM 26 


-5.0 ~ 1.0 


29.9 


0.089 


0.022 


0.0065 


0.29 


7.5 x 10~ 8 


2.2 x 10~ 6 


0.73 


W 


EL 29 


-5.0 ~ 1.0 


21.6 


0.083 


0.11 


0.025 


0.38 


3.0 x 10~ 7 


6.5 x 10~ 6 


1.3 


N 


Red Lobe 


VLA 1623 


6.5 ~ 13.0 


22.7 


0.096 


0.25 


0.056 


0.41 


6.0 x 10 -7 


1.4 x 10 -5 


2.9 


N-E 


EL 32 


6.5 ~ 10.0 


5.7 


0.17 


0.50 


0.028 


2.9 


1.7 x 10~ 7 


9.7 x 10~ 7 


0.049 


S-E 


LFAM 26 


6.5 ~ 12.0 


19.4 


0.11 


0.25 


0.048 


0.55 


4.5 x 10~ 7 


8.6 x 10~ 6 


1.5 


E 


LFAM 26 


7.0 ~ 12.0 


24.2 


0.0089 


0.037 


0.089 


0.36 


1.0 x 10~ 7 


2.5 x lO" 6 


0.55 


W 


EL 29 


6.5 ~ 12.0 


21.6 


0.13 


0.22 


0.047 


0.60 


3.7 x 10~ 7 


8.0 x 10~ 6 


1.5 


s 



Note. — The CO ( J = 3 — 2) emission lines are assumed to be optically thin. Therefore, the physical quantities listed here give the lower limits. 

a Relative Position of the lobe from the driving source: N, S, E, and W mean north, south, east, and west, respectively. For example, N-E denotes the lobe 
flowing from the driving source toward north-east direction. 



Table 3. Molecular Outflow Parameters Obtained From CO ( J = 1 - 0) 



Name 


Velocity Range 
(km s -1 ) 


Vout 

(km s — 1 


Rout 
) (PC) 


M out 
(1O- 2 M ) 


Pout 

(M km s" 1 ) 


^dyn 

(ioV) 


Mout 

(M yr- 


L ) 


Pout 

(M km s~ 1 yr _ 


Pout 

i) (io- 2 p ) 


Position 


Blue Lobe 


VLA 1623 


-3.0 


~ 1.0 


20.8 


0.52 


3.9 


0.78 


2.5 


1.5 x 10" 


-6 


3.2 x 10- 


-5 


5.8 


S-E 


VLA 1623 


-3.0 


~ 1.0 


18.8 


0.084 


0.34 


0.065 


0.44 


8.0 x 10" 


-7 


1.5 x 10- 


-5 


2.3 


N-W 


EL 32 


-3.0 


~ 0.5 


9.4 


0.22 


10.0 


0.94 


2.3 


4.4 x 10" 


-6 


4.1 x 10- 


-5 


3.3 


N-W 


LFAM 26 


-3.0 


~ 1.5 


20.7 


0.053 


0.023 


0.0050 


0.25 


9.5 x 10" 


-8 


1.9 x 10- 


-6 


0.33 


E 


EL 29 


-3.0 


~ 1.5 


13.5 


0.074 


0.25 


0.033 


0.54 


4.7 x 10" 


-7 


6.2 x 10- 


-6 


0.69 


N 


IRS 44 


-3.0 


~ 1.0 


8.7 


0.036 


0.15 


0.013 


0.40 


3.7 x 10" 


-7 


3.2 x 10" 


-6 


0.26 


S-W 


Red Lobe 


EL 32 


6.5 ~ 


9.0 


4.6 


0.31 


3.2 


0.15 


6.5 


5.0 x 10" 


-7 


2.3 x 10- 


-6 


0.087 


S-E 


LFAM 26 


6.5 ~ 


9.0 


15.8 


0.10 


0.31 


0.048 


0.64 


4.8 x 10- 


-7 


7.6 x 10- 


-6 


0.99 


E 


LFAM 26 


6.5 ~ 


9.0 


15.4 


0.068 


0.045 


0.0070 


0.43 


1.0 x 10- 


-7 


1.7 x 10- 


-6 


0.21 


W 


EL 29 


6.5 ~ 


9.0 


17.6 


0.096 


0.34 


0.062 


0.54 


6.5 x 10- 


-7 


9.7 x 10" 


-6 


1.7 


S 


EL 29 


6.5 ~ 


9.0 


16.9 


0.083 


0.21 


0.034 


0.48 


4.3 x 10" 


-7 


7.2 x 10- 


-6 


1.0 


N 


IRS 44 


6.5 ~ 


9.0 


8.4 


0.089 


0.14 


0.012 


1.0 


1.3 x 10" 


-7 


1.1 x 10" 


-6 


0.083 


N 



Note. - 



The CO ( J = 1 — 0) emission lines are assumed to be optically thin. Therefore, the physical quantities listed here give the lower limits. 
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Table 4. Molecular Outflow Parameters Compiled From the Literature in the p Ophiuchi Main 

Cloud 



Name 


Vout 


Afout 




^out 


Reference a 




(km s" 1 ) 


(M e ) 


(Mq km s -1 yr _1 ) 


(LQ) 




WL 6 


7.5 


1.5 x icr 3 


2.0 x 10~ 6 


1.2 x 10 -3 


1 


WL 10 


7.5 


5.7 x 10" 3 


3.5 x 10" 6 


2.2 x 10" 3 


1 


AN 


3 


8.2 x icr 3 


2.6 x 10~ 6 


7.0 x 10~ 3 


2 


AS 


2 


2.9 x icr 3 


3.7 x 10~ 6 


2.9 x 10~ 3 


2 


MMS126 


(5.5 


1.0 x 10" 3 


0.19 x 10" 6 


0.059 x 10~ 3 


3 


WL 12 






1.4 x 10~ 6 




4 


IRS 43 






3.4 x 10~ 6 




4 


IRS 48 






6.9 x 10~ 6 




4 


IRS 51 






1.3 x 10~ 6 




1 



a Reference: 1. ISekimoto et al ] l|l997l l. 2. | Kamaz aki et al ] &200$ ). 3. IStanke et al] l|2006h . 4. iBontemps et alJ l|l996h 

Note. — All the quantities listed in the table arc corrected to D = 125 pc under the assumption that observed lines 
are optically thin. The values of outflow mechanical luminosities for WL 5, and WL 10 are estimated from the expression 

Lout = -FoutVout/2. 



Table 5. Global Properties of Nearby Cluster Forming Clumps 



Name 


Distance 


Mass 


Radius 


Velocity Width Virial Ratio 


^turb 


Ltot -^grav 


-Ftot 




(pc) 


(Mq) 


(PC) 


(km s" 1 ) 


(£©) 


(Mq km s — 1 yr ) 





Serpens 1 260 210 0.46 2.0 1.1 0.14 - 0.32 1.3 1.9 X 10" 4 5.2xl0~ 4 

pOph 2 125 883 0.8 1.5 0.22 0.06 - 0.12 0.2 1.3 x 10" 3 1.2 x 10~ 4 

NGC 1333 3 220 573 0.44 2.8 0.76 0.41 - 0.90 0.65 1.8 x 10" 3 6.0 x 10~ 4 



Note. — Reference: 1 Sugitani ct al. (2010]), 2 this paper, ^Knee & : Sandcll (2000). For all the clumps, the contribution of magn etic 
support is taken into account for estimation of F gIa , v whose value is reduced by a factor of about 2. See ISugitani et~ai1 feoich m 
detail. The outflow parameters listed in the table are the lower limits because they are estimated under the assumption that the lines 
are optically thin. 
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Kkm/s 




RA (J2000) 

Fig. 1. — CO (J = 3 — 2) total integrated intensity map in the velocity range from vlsr = —9.8 km 
s _1 to +12.2 km s" 1 toward the p Ophiuchi main cloud. The CO ( J = 3 — 2) integrated intensity 
is shown in color in units of K km s _1 . The circles indicate the positions of B stars SI, DoAr 21, 
and SR3. 
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Kkm/s 




28rn0s 27m0s 16h26m0s 

RA (J2000) 



Fig. 2. — CO ( J = 1 — 0) total integrated intensity map in the velocity range from vlsr = —3.2 
km s _1 to +8.0 km s _1 toward the p Ophiuchi main cloud. The CO ( J = 1 — 0) integrated intensity 
is shown in color in units of K km s^ 1 . The circles indicate the positions of B stars SI, SR3, DoAr 
21, and HD 147889. 
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Fig. 3. — (a) CO (J = 3 — 2) map of the p Ophiuchi main cloud on the 850 ^m image obtained by 
Johnstone et al.l (|2000i ). Blue contours represent blueshifted CO ( J = 3 — 2) intensity integrated 



between —9.8 km s _1 and 2.2 km s , starting from 18 K km s _1 at intervals of 4 K km s . Red 
contours represent redshifted CO ( J = 3 — 2) intensity integrated between 7.4 km s _1 and 11.8 km 
s , starting from 1.6 K km s _1 at intervals of 0.8 K km s _1 . The gray scale shows the 850 /j,ui 
image in units of Jy/beam. The dense subclumps are designated by A, bl, B2, C, E, and F. (b) 
CO ( J = 3 — 2 ) map of the p Op h iuchi main cloud overlaid with the positions of the embedded 
YSOs listed by Ivan Kempen et al.l (|2009i ) and some other YSOs with the yellow circles. Positions 
of several peaks of outflow components detected by the CO ( J = 3 — 2) observations are indicated 
by the alphabets (a) through (j). 
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Fig. 4. — CO (J = 1 — 0) map of the p Ophiuchi main cloud on the 850 image obtained by 
Johnstone et al.1 (12000 ) . Blue contours show the CO ( J = 1 — 0) intensity integrated between —4.0 
km s _1 and 1.2 km s _1 , starting from 8.9 K km s" 1 at intervals of4.1 K km s _1 . Red contours 
show the CO ( J = 1 — 0) intensity integrated between 6.4 km s _1 and 9.2 km s _1 , starting from 
4.1 K km s _1 at intervals of 2.8 K km s _1 . The gray scale, symbols, and alphabets are the same 
as those of Figure [3l 
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Fig. 5. — CO (J = 3 — 2) spectra toward several high-velocity lobes indicated by alphabets in 
Figure Each spectra is smoothed in a 1' x 1' grid centered at the position of the peak intensity. 
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Fig. 6. — CO (J = 1 — 0) spectra toward several high-velocity lobes indicated by alphabets in 
Figure [H Each spectra is smoothed in a 1' x 1' grid centered at the position of the peak intensity. 
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Fig. 7. — CO ( J = 3 — 2) map of the p Ophiuchi main cloud. The contours are the sa me as those of 
Figure [3l The crosses and plu s es ind icate the positions of the H2 knots identified by iGomez et al. 
(|2003l ) and iKhanzadvan et al.l ()2004l ). respectively. 
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Fig. 8. — (a) CO (J = 3 — 2) map of the p Ophiuchi main cloud overlaid on the H 13 CO + (J = 1 — 0) 
integrated intensity map obtained by lMaruta et al.1 (120101 ) in gray scale, (b) C O (J = 3 — 2) map o f 
the p Ophiuchi main cloud. The positions of the H 13 CO + cores identified by lMaruta et al.1 (201CJ) 
are overlaid by the squares. For both the panels (a) and (b), the blue and red contours are the 
same as those of Fig. [3l 
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Fig. 9. — Line-width-radius relation of the H 13 CO + cores identified bv lMaruta et al.l ( 201(3 ) . Filled 
circles are the cores that appear to be associated with the EL 32 outflow lobes, while open circles 
the other cores. 



